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Abstract: We report on the detection of HCO^ and ^'^CO emission in the rotational transition J=l-0 
in the vicinity of the shock front at the southern border of the supernova remnant RCW 103, where 
previous infrared observations suggest an interaction with a molecular cloud. The observations were 
carried out with the Australian Millimeter Radiotelescope at Mopra. We observed a depletion of HCO"'' 
behind the supernova shock front. In addition, we studied the interstellar medium over an extended 
region towards RCW 103 based on archival A 21 cm HI line observations from the Australia Telescope 
Compact Array (ATCA) and the Parkes Telescope. No atomic gas is observed in emission in coincidence 
with the molecular feature. This absence is interpreted in terms of self absorption processes. 
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1 Introduction 

Supernova remnants (SNRs) mark the catastrophic death 
of massive stars, and they modify irreversibly the sur- 
rounding medium into which the shock front is ex- 
panding. Roughly half of all G alactic SNRs are lo- 
cated near large molecular clouds l)Huang fc ThaddeusI 
1198611 ■ presumably the parent clouds of the massive 
stellar progenitors. When an SNR interacts with a 
molecular cloud, the shock wave depletes and excites 
molecular species, significantly altering the local chem- 
istry. 

RCW 103 (C332.4-0.4) is a Galactic SNR that shows 
clear evidence of an interaction with a molecular cloud 
near its southern limb. At radio wavelengths, this SNR 
appears as an almost complete cir cular shell with a di- 
ameter of 8' llCaswell et al l llQSOH . Optical filaments 
are s een towards the brightest regions of t he radio 
sheU (Ivan den Bergh et al.llT97l lRuizlll983D. Based 
on HI absorption measurements, ICaswell et al.l lll975f) 
estimated a d istance of 3.3 kpc f or RCW 103. Using 
this distance, ICarter et al.l il997ri studied the expan- 
sion of the optical filaments and calculated an age for 
the remnant of about 2000 years. In X-rays, RCW 103 
presents a structure compatible with shell S NRs in the 
youn g double-shock stage of their evolution l)Chevalieil 

Figure 1 shows a radio continuum imag e of RCW 
103 at 1.4 GHz from iR.evnoso et alJ lj2004ri (left) and 
a soft X-ray image (right) extracted from the Chandra 
Supernova Remnant Catalog. Good concordance is ob- 
served between the non-thermal radio emission and the 
soft X-ray emission. A compact central object (CCO) 
has been detected in hard X-rays JGotthelf et alll99gll 
near R.A. 16''17™35fi, dec. -51° 02' 25" (J2000) with 
no detected radio counterpart. There are strong indi- 



cations that this CCO is pa rt of a binary system (e.g., 
IBecker fc Aschenbachll2002t) . 

Observations in the 2. 122 ^tm line of the H2 and 
other near infrared lines jBurton fc Spvromiliclfl993l : 
lOliva et alJll99(ll . suggest that the remnant is inter- 
acting with molecular material on the southern side. 
Burton and Spyromilio conclude that the H2 spectrum 
of RCW 103 is typical of dense gas (> 10* cm'^) col- 
lisionally excited by low v elocity (10 - 20 km s~^ ) 
shocks. In a recent paper, IChevalie3 1)200511 suggests 
that RCW 103 is not embedded in a molecular cloud, 
but rather that the blast wave is interacting with previ- 
ous circumstellar material swept up by the progenitor 
star. 

In a search for shocked molecular gas near the 
southern border of RCW 103 we carried out observa- 
tions in two molecular lines: ^^CO J=l-0 and HCO'*' 
J=l-0. The first one is the brightest molecular line, 
very sensitive to the presence of molecular gas. Its de- 
tection confirms the existence of molecular concentra- 
tions on the southern side of RCW 103. HCO"*" emis- 
sion is enhanced at high densities and this molecule 
is particularly sensitive to shock emission, providing a 
good tracer of the shock location. However, this line 
is generally quite weak. 

The interstellar HCO"*" ion has been extensively 
studied in the past. In well-shielded regions, where 
interstellar radiation plays a minor role (molecular re- 
gions), the dominant reactions that produce HCO"'' 
are the following: H2 -I- Cosmic Rays Yi^ + H + 
e", then H+ Ha ^ H3+, and finally H3+ + CO 

HCO^ + H2. There is much controversy in the 
literature about t he way in which shock s modify the 
HCO+ chemistry. Ilglesias fc Sil^ il978D showed that 
HCO'*' must be depleted behind slow shocks. However, 
HCO^ observations towards IC443 showed an abun- 
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Figure 1: Left: Non-thermal radio emission from 
RCW 103 at 1.4 GHz (l^eynoso et al. 2004). The 
contour levels are 0.1, 0.2, 0.3, 0.4 and 0.5 Jy 
beam""'^. Right: Soft X-ray emission in the energy 
band 1.6-3.0 keV {Chandra Supernova Remnant 
Catalog) 
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dance enhancement of this molecule behind t h e super- 
nova shock front (e.g.. lDickinson et alj|l980l) . lElitzuil 
ljl983l ^ proposed that this enhancement resulted from 
increased ionization due to cosmic rays trapped by 
shock-compressed magnetic fields. However, subse- 
quen t observations of IC443 ijZiurvs et alJl98fll : lvan Dishoeck et alJ 
Il993f) showed that the HCO^ abundance in the shocked 
gas is not significantly enhanced compared with pre- 
shocked values. 

In this paper we investigate the interplay between 
RCW 103 and its surrounding medium based on obser- 
vations of ^^CO and HCO'*'. A complementary analy- 
sis of HI A 21 cm data is also presented. 



Figure 2: Left: Radio continuum emission at 
1.4 GHz in an exten ded field around RCW 103 
(jRevnoso et al.ll2004|) . The nine dots towards the 
South of the SNR indicate the pointings for the 
molecular line observations. The grayscale is in 
units of Jy beam~^. The beam, 50" x 50" , is 
plotted in the bottom right corner. The rms noise 
level of this image is 5.5 mJy beam~^. Right: A 
line profile at constant R.A. (16''17'°37fi) for 1.4 
GHz radio continuum and X rays (1.6-3.0 keV, 
from Chandra Supernova Remnant Catalog). The 
crosses with numbers on top indicate the pointings 
of the molecular line observations. 



2 Observations 



the shell centre (smaller negative declinations) on the 
left side of the plot. The crosses with numbers at the 
top indicate the positions of the molecular line obser- 
vations, and correspond to the numbered profiles in 
Figure 3. Pointings 1 and 2 lie inside the SNR, point- 
ing 3 is at the rim (as seen in radio wavelengths), and 
the rest probe gas external to the remnant. 



2.1 Millimetre Line Observations 

We observed the J = 1 - rotational transitions of ^^CO 
{ly = 115 GHz) and HCO+ {ly = 89 GHz) using the 22- 
m Mopra antenna, near Coonabarabran, NSW, Aus- 
tralia. Spectra were recorded in position-switching 
mode, with the telescope alternating between the tar- 
get and an emission-free reference position. The back- 
end consisted of an auto-correlator, configured to have 
a bandwidth of 64 MHz divided into 1024 channels, 
which provided a velocity resolution of ~0.2kms^^ 
over a usable bandwidth of ~ 120 kms^^. The pointing 
accuracy was checked before each mapping run by ob- 
serving a nearby SiO maser, and is estimated to be bet- 
ter than 10". Nine pointings, at constant right ascen- 
sion (R.A. — 16'^17™37fi) were observed towards the 
southern rim of RCW 103 with a resolution of ~34". 
Pointings were spaced at 30" intervals in declination, 
from -51°05'30" to -51°09'30" (Figure 2, left). Also 
in Figure 2 (right), we include a radio continuum and 
an X-ray emission profile at constant Right Ascension 
(R.A. = 16''17™37fi) to show the SNR shock front po- 
sition. In this picture, the SNR outer shell peaks at 
dec. ~ —51° 6', and the intensity decreases towards 



2.2 HI Observations 

The A 21 cm HI da ta used in the presen t paper are 
described in detail in lR.evnoso et all ||2q3). These ob- 
servations were carried out using the Australia Tele- 
scope Compact Array (ATCA) in the 750A array of a 
1° X 1° region around RCW 103, during a 12 h run 
on 2002 January 22. The baselines in this array range 
from 76.5 to 735 m. A second observation of 12 h with 
the telescope in the EW 367 array (baselines from 46 
to 367 m) took place on 2002 April 1. The angular 
resolution of the combined data is 50". A correlator 
configuration of 1024 channels was used with a total 
bandwidth of 4 MHz centered at 1420 MHz. The ve- 
locity resolution at this frequency is 1 km s^^. Contin- 
uum data were obtained simultaneously with a band- 
width of 128 MHz centered at 1384 MHz. Single dish 
data from the Southern Galact ic Plane Survey (SGPS) 
(iMcClure- Griffiths et al.]|2005D in the A 21 cm HI line 
were added in order to recover the contribution from 
lower spatial frequencies. 
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3 Results 

3.1 Molecular Gas 

Figure 3 displays all the observed profiles in the J=l- 
transition of ^'^CO and HCO^ in the velocity range 
between —60 km s~^ and —20 km s~^. In this ve- 
locity interval, the ^^CO spectra have several com- 
ponents, with the brightest one centered near —48 
km s~^. The HCO"^ spectra are, as expected, sim- 
pler with only one conspicuous feature centered at a 
velocity varying between —45 and —50 km s~^ along 
the pointings, broadly in coincidence with the most 
intense '^^CO peak (all kinematical velocities in this 
paper are referred to the LSR). 

In Table 1 we present the observed parameters 
of the molecular line components peaking near —48 
km s~^. The peak temperatures (Tpeak), central ve- 
locities (Vc) and velocity widths (Av) were obtained 
from Gaussian fits to the main components of the spec- 
tra. The ^^CO central velocities correspond to the 
most intense peaks. The rms noise is ctco ^ 0.15 K 
and crjjco+ ~ 0-02 K. The last column represents the 
ratio between the HCO"^ and ^^CO peak temperatures. 

3.2 Atomic Gas 

Figure 4 shows the distribution of HI in grayscale from 
V — —75 km s~^ to V = +i3 km s~^ integrated in 
10 km s~^ steps. The contours correspond to the ra- 
dio continuum emission of RCW 103. These velocity 
slices reveal the patchy appearance of the cold gas in 
the surroundings of RCW 103. The hole in the emis- 
sion evident in the map at —21.3 km s~^, corresponds 
to absorption by the Sagitarius-Carina Galactic arm, 
whose distance is ~1.5 kpc. Since HCO"'" is observed 
only between —52 and —43 km s~^. Figure 5 shows 
the HI A 21 cm emission integrated over this velocity 
range to determine if there is an atomic gas counter- 
part to the molecular feature. There is a hint of a 
shell of atomic gas open to the NE of the SNR, and 
a structure partially overlaping the remnant towards 
the East. Assuming that the gas is optically thin and 
using the relation N(HI) = 1.82 x 10^** / Tb dv, where 
Ts is the brightness temperature, the HI column den- 
sity integrated in the above velocity range is estimated 
to be N{HI) ~ 1.6 x 10^^ cm"^ for the shell and the 
East clump. Based on these values and assuming a dis- 
tance of 3.3 kpc for RCW 103 (from section 4.2), we 
calculate an HI mass --^ 800 Mq for the East clump 
and aproximately 4400 Mq for the open shell. Curi- 
ously, no atomic gas is detected close to the southern 
border of RCW 103 down to a detection level of 3a". 
This fact is discussed later. 

4 Discussion 

4.1 Molecular Gas 

From Figure 3, it is clear that the HCO^ intensity de- 
creases from outside to the interior of the SNR shell, 
with the emission being weakest for the pointings 1 to 
3. The coincidence between the shock front position 



and the change in the HCO intensity strongly sug- 
gests the existence of a physical interaction between 
the SNR and the molecular gas. To investigate if 
such a change is the consequence of the passage of the 
shock front or is simply due to a decrease in the am- 
bient molecular density, it is necessary to compare the 
HCO"^ emission with that of another molecule that 
remains unnaffected by the shock passage. For that 
purpose, we compared the HCO^ emssion with that of 
^^CO. The use of ^^CO as a reference molecule may be 
problematic because the excitation conditions between 
'^^CO and a high dipole moment species like HCO"'' 
could be different. However, assuming that ^^CO and 
HCO"^ are uniformly mixed, the comparison would be 
meaningful. We find that the t^'^^^/t'^'^^ ratio is 
lower for pointings 1 to 3 than for those that lie well be- 
yond the SNR shock front. For the inner pointings, the 
lower values in the intensity of the HCO^ emission are 
probably indicative of shock depletion of this molecule. 
Thus, our results confirm that the SNR shock affects 
the chemistry of the molecular gas, but do not sup- 
port models that predict HCO"'' enhancement by the 
passage of a shock. 

However, the decrement is not as low as that pre- 
dicted by the Iglesias & Silk (1978) model, which pro- 
duces a decrease of two orders of magnitude in the 
HCO"'' abundance behind a shock front, on time scales 
of 1000 yr and for shock velocities of ~ 10 km s^^. A 
short interaction time or a shock front with different 
velocity could explain the smaller decrease observed in 
the HCO"^ emission behind the shock, as compared to 
the theoretical predictions. Nevertheless, to determine 
the HCO"'' abundance accurately, it would be useful 
to observe a higher dipole moment m olecule, such as 
HCN, rather than ^^CO (see lZiurvs eTal-lOSftl for the 
reference molecule. 

4.2 Kinematical Distance 

Previous distance studies suggested a range of pos- 
sible values between 3.1 and 4.6 kpc for RCW 103 
iRevnoso et all2004f) . The velocity at which the HCO+ 
emission is produced, coincident with that of the most 
intense ^^CO peak, allows us to constrain the systemic 
velocity of RCW 103 to about -48 km s~\ This ve- 
locity corresponds to a dist ance of ~ 3 . 3 kpc using the 
Galactic rotation model of iFich et al.l (|1989D . 

4.3 Atomic Gas 

From Figure 4 and Figure 5, it is clear that on the 
southern side of RCW 103, where the remnant may 
be interacting with a molecular cloud, we do not ob- 
serve dense HI emission adjacent to the SN shock. It 
is possible that the lack of HI emission in this region is 
not due to the absence of atomic hydrogen, but to self- 
absorption processes, as can be seen in Tycho's S NR. 
Proper motion measurements (iRevnoso et alll997f) showed 
that the lowest shock expansion velocities around the 
outer shell of Tycho's SNR are found to the East, and 
it was suggested that a higher ambient density was de- 
celerating the exp ansion in this directi on. However, a 
subsequent study iRevnoso et al.iri999ri detected faint 
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Figure 3: The ^^CO and HCO+ spectra are displayed in pairs with the ^^CO profiles in the upper panel 
and HCO+ in the lower one for every pointing, numbered as in Figure 2. The plotted velocity range is 
[—60 km s^^, —20 km s~^] for both molecular transitions. The temperature scales vary from — 1 to 8 K 
for i^CO and from -0.1 to 0.4 K for HCO+. 



Table 1: Observed and derived parameters. 
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Figure 4: HI images between —75 and +43 km s^^, each one integrated over ^ 10 km s^^. The velocity 
shown in the upper left corner of each image corresponds to the first channel of integration. The grayscale 
ranges from 400 to 1200 K km s~^. A few contours of the radio continuum of RCW 103 are included for 
comparison. 



HI emission in this region. The observations can be 
explained if a cooler HI comp onent is self-absorbing a 
warmer one. iLee et aJ] ||2Q^) confirmed this hypoth- 
esis with detection of the molecular emission counter- 
part of the HI self- absorption. The presence of am- 
bient, unshocked molecular gas is indicative of low 
temperatures. A similar scenario for RCW 103 also is 
consistent with Burton & Spyromilio's (1993) model, 
in which a reverse shock is propagating back into the 
swept-up atomic gas following the interaction of the 
blast wave with a dense molecular cloud. The blast 
wave is propagating at 15 km s~^ and the molecular 
cloud is located ahead of the neutral hydrogen mass. 
Moreover, since the reverse shock has already overrun 
the HI cloud, the shocked gas will be much warmer 
than the unshocked gas which is mixed with the molec- 
ular ambient cloud, thus producing the observed self- 
absorption. In this context, the eastern HI clump 



could be part of an extended structure which has a 
uniform temperature and for which the emission is not 
self-absorbed. A sketch of the possible gas distribution 
is shown in Figure 6. The cool molecular component 
appears closer to us than the warmer atomic gas and 
the supernova shock front is currently running into the 
molec ular cloud. 

In lRevnoso et all (l2004fl . the absorption spectrum 
towards RCW 103 was analyzed by applying different 
filtering techniques. In their Figure 2, an absorption 
feature that was seen well above the noise (> 6a) at 
-(-34 km s^^, disappeared after deeper filtering and 
hence, was considered not to be caused by absorption 
against the continuum emission from RCW 103. To 
review a possible association with RCW 103, a map of 
the HI distribution at this velocity was produced. 

Figure 7 shows the HI distribution integrated within 
4 km s~^ around -f34 km s~^. The neutral gas depicts 
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Figure 5: Total HI A 21 cm emission integrated between —52 and —43 km s ^. The grayscale varies 
between 650 and 1100 K km s-\ Contours at 800, 900, 1000 and 1050 K km s^^ are included to 
emphasize the morphology of the atomic gas. A few contours from the radio continuum image of RCW 
103 arc shown for comparison. The rms noise level of the image is 47 K km s^^ and the resolution is 
50" X 50". 
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Figure 6: A schematic model of the RCW 103 in- 
teraction with a molecular cloud to the South. 



an open shell-like structure well-matched to the SNR 
ring. Since the molecular data presented here and the 
HI absorption study in Reynoso et al. (2004) have 
shown that the systemic velocity of RCW 103 is —48 
km s~^, it is highly unlikely that this HI shell at -1-34 
km s~^ (which corresponds to a kinematic distance of 
19 kpc using the Galactic rotation model of Fich et al. 
1989) has been driven by the SN shock front. 

Although there is a good morphological match be- 
tween the radio continuum extent and the HI feature, 
it is very difficult to explain any relationship between 
the SNR and this HI structure. The most likely expla- 
nation is chance juxtaposition and no relationship. 




1620 19 18 17 16 15 
RIGHT ASCENSION (J2000) 



Figure 7: HI emission integrated within 4 km 
around -1-34 km s~^. The grayscale varies between 
30 and 70 K km s^^. The HI contour levels (in 
black) are 30, 40, 50, 60 and 65 K km s"^ A few 
contours (in grey) representing the radio contin- 
uum of RCW 103 arc shown for comparison. 
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5 Summary 

We have investigated the interstellar medium around 
the SNR RCW 103. The present ^^CO J=l-0 and 
HCO"'' J=l-0 observations provide new evidence for 
an interaction between the SNR and a molecular cloud. 
We also investigated the surrounding atomic gas based 
on the A 21 cm HI line. The main results can be sum- 
marized as follows: 

(a) From J=l-0 ^^CO and HCO+ observations to- 
wards the southern side of RCW 103, we find that the 
HCO"'' spectra have only a single component which has 
at the same kinematical velocity of ~ —48 km s~^ as 
the brightest feature observed in the '^^CO spectrum. 

(b) The T"CO+^rpi2co ^j^j^ velocity range 
has lower values for the interior pointings (1 to 3) than 
for those that are well away from the SNR shock front. 
The low ratio values for the innner positions may be 
due to HCO^ depletion caused by the remnant shock 
front colliding with the molecular gas. 

(c) The molecular observations allows us to con- 
strain the systemic velocity of RCW 103 to about —48 
km s~^. This velocity corresponds to a distance of ~ 
3.3 kpc for the re mnant, using the Galactic rotation 
model of lFich et al.. (.IQSa'l . 

(d) We do not detect HI emission on the southern 
side of the remnant at —48 km s~^, where molecular 
gas is observed. We suggest that this result is due to 
self-absorption processes, similar to those observed in 
Tycho's SNR. 
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